Polymer chain in a quenched random medium: 
slow dynamics and ergodicity breaking 
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O \ The Langevin dynamics of a self - interacting chain embedded in a quenched random medium is 

investigated by making use of the generating functional method and one - loop (Hartree) approx- 
^"*^ ■ imation. We have shown how this intrinsic disorder causes different dynamical regimes. Namely, 

within the Rouse characteristic time interval the anomalous diffusion shows up. The corresponding 
subdiffusional dynamical exponents have been explicitly calculated and thoroughly discussed. For 
the larger time interval the disorder drives the center of mass of the chain to a trap or frozen state 
provided that the Harris parameter, (A/b d )N 2 ~ vd > 1, where A is a disorder strength, b is a Kuh- 
nian segment length, N is a chain length and v is the Flory exponent. We have derived the general 
equation for the non - ergodicity function f(p) which characterizes the amplitude of frozen Rouse 
modes with an index p — 2nj/N. The numerical solution of this equation has been implemented and 
shown that the different Rouse modes freeze up at the same critical disorder strength A c ~ TV -7 
where the exponent 7 » 0.25 and does not depend from the solvent quality. 



PACS numbers: 61.25.Hq,78.55.Qr,66.90.+r 
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I . I. INTRODUCTION 

a 

The statical and dynamical properties of polymer chains in a random medium bares still many open relevant 
questions in the current discussion on polymer physics and disordered systems [Q. This is due to the fact that this 
model naturally occurs under the interpretation of many experiments: polymer adsorption on heterogeneous surfaces 
' pL polymer diffusion in the swollen polyvinylmethyiate gels ||, electrophoresis ||), etc. 
. The theoretical investigation of the problem has been stimulated by the Monte - Carlo simulation on a site - 
diluted lattice Q. These results essentially extended to the dynamic properties (anomalous diffusion, chain length, 
0^ \ N, dependence of the diffusion coefficient D, etc.) in the number of papers || |7|, |[ [|. The authors of these papers 
nI 1 have found new dynamic laws where chain diffusion is slower than even reptation. Such behavior emerges from 
the presence of "entropic traps" , i.e. regions which are relatively free from obstacles and hence entropically more 
preferable. The diffusion is slowed down substantially by the presence of narrow channels between traps, so that the 
chain is forced to squeeze through them. 

A large number of theoretical studies emphasizing the statical [if], [f^, |f3|, Q |R| [l6], [l7] [L8| [ll| and dynamical 
fi"2| . po| , pi| , |2^ , |23| , |24| , |2E| , p6| aspects of the problem. When discussing the static properties one should discriminate 
between the annealed and quenched averaging. If the system is large enough and the chain is still rather mobile, 
so that it experiences different disorder environment (during the time of experiment) the problem is reduced to the 
annealed one. In this case the only effect of disorder (which is characterized by the second moment A of the quenched 
^ ■ potential fluctuation) is the reduction of the excluded volume parameter v, i.e. v — > v — A ||,||||,||. In the mean 
time it is well known [|l4| |lj, |2(], ^2|, that the renormalization group (RG) consideration for the discussed problem 
suffers from he lack of a stable fixed point for the renormalized second moment Ar, so that the perturbation theory 
can not be used. Nevertheless, it was argued 22, 2^, Q that for the weak disorder in the thermodynamical limit the 



quenched and annealed averaging lead to the same results. The authors of ref. |f^, [ijj drew similar conclusions by 
studying a self - avoiding chain on dilute lattices and making use the polymer - magnetic analogy and n — * trick j27[ . 
On the other hand it has been shown in ref. |f2|, [ll| that "entropic traps" caused by disorder affect the equilibrium 
spatial distribution Zn of the polymer chain strongly for d < 4. In this case the distribution is inhomogeneous, so 
that typical and average values of Zn are different. Under these conditions the motion of the chain center of mass 
can be trapped. The corresponding diffusion coefficient, D, is scaled with the length of chain N as D ~ exp(— AN a ), 
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where the "specific heat" exponent a — 2 — vd and v is the Flory exponent. It can be seen that for rather long chains 
(remember that a > at d < 4) the diffusion is effectively suppressed and the whole system becomes nonergodic. 

This is an indication of the multiple local minima configurations which could manifest itself through the replica 
symmetry breaking (RSB) scheme |2£| |3(J. It have been shown there that within the RG - approach the conven- 
tional replica symmetric (RS) fixed point (which used to be considered as providing new universal disorder - induced 
critical exponents) are unstable with respect to an RSB - solution. In this case the structure of the renormalized 
interactions develops strong RSB and the values of the corresponding interaction parameters are getting large. 

In this paper we use an alternative, dynamical, approach in order to treat nonergodic regimes for the chain in a 
quenched random medium. It is of interest that the nonergodic, glassy - like regimes have been recently seen in quasi 
- elastic neutron scattering experiments on dynamic properties of flexible polymer chains filled with nano - particles 
(hydrophilic aerosil) |3lJ . The authors of |3l| have seen the enhancement of elastic component of neutron scattering 
from chains whose dynamics is getting strongly restricted and slowed down upon the growth of nano - particles 
fraction. To the best of our knowledge the theoretical studies of these regimes is lacking. We will use the Langevin 
dynamics and Martin - Siggia - Rose generating functional method as well as the Hartree approximation in order 
to derive the equation of motion for our random model. We will show that at time t — > oo this leads to the Gotze - 
like equation B3| for the non - ergodicity function, f(p), which describes frozen Rouse - mode p - dependent states. 
The dynamic phase diagram in terms of Rouse mode indices and disorder parameter A is numerically calculated and 
thoroughly discussed. 



II. PRELIMINARIES 



A. Model 



The chain conformation is characterized by the d-dimensional vector - function R(s, t) of the time t and s (where 
1 < s < N), which labels beads of the chain. Besides intrachain interactions the chain experiences a quenched random 
external field, l^{R(s, t)}, so that the whole Hamiltonian has the form 



H = \ e Y. [V S R(M)] 2 + #i„t {R(M)} + ^{R(M)} , (2.1) 

1 3=0 

where e = dksT jb 2 is the elastic modulus with the Kuhn segment length b, N is the length of chain and the finite 
difference V s Rj(s, t) = Hj(s + !,<) — Rj(s, i) . The intra - chain Hamiltonian reads 



H iQt {R(s,t)} = i]T £>(R(M)-R(a',t)) 

s=0 s'=0 
N-1N-1N-1 



s=0 s'=0 s"=0 



where v(r) and u>(ri,r2) are the second and third virial coefficients correspondingly. The quenched random potential 
V{R(s,t)} is assumed to be Gaussian distributed, i.e. 

(V(r)V(r')) = A^(r-r'), (2.3) 

where dispersion A is one of the main control parameters of the problem and has dimensions of volume (in these units 
of measurement ksT = 1). This will become important below. Eqs.(2.1) - (2.3) represent and specify our model. 



B. Annealed and quenched disorder 



Let us consider first the case when the characteristic times for chain's configurations and chain's center of mass 
positions in a disordered medium are of the same order of magnitude. Then in the course of an experiment the chain 
experiences all possible quenched field realizations. This corresponds to the annealed disorder and the corresponding 
free energy i^nneai = — hi (S)y, where S is the partition function at a given realization of V(r) and (. . . ) v is the 
averaging over the field V(r) distribution. It can be shown that in this case 



<S) y = So{«-A} 



(2.4) 
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where 2o{^} is the partition function of the pure (i.e. without disorder) system with the second virial coefficient 
v. As a result the only effect of the disorder is the reduction of the second virial coefficient, i.e., a reduction of the 
excluded volume. However, it may change the sign of the effective second virial coefficient, and can cause therefore 
collapsed states of the chain. 

Nevertheless, in a medium with a strong disorder the chain is preferably trapped in some regions where the depth 
of the quenched random potential exceeds ksT. In this case the chain is pinned down in some particular place of a 
disordered medium and experiences only a local quenched field. This corresponds to the quenched disorder and the 
relevant free energy -Fq Uenc h = — (lnS)y. For the site - diluted lattice medium model it was arg ued by Machta @ [f| 
that while the size of the chain is unaffected by the disorder (i.e. R ~ bN v ) the whole spatial distribution of the 
chain is correlated with the disorder for d < 4. The effect of the disorder shows up as an essential singularity in the 
N dependence of the so - called typical value of the partition function, S typ = exp (— F qucnc h), as well as the chain's 
center of mass diffusion coefficient D. 

The chain with a frozen center of mass explores only a local region of the medium. As a result a simple estimation 



of the typical value of the partition function reads [|35| 



"typ 



-oMexp^YV--/ 2 } , (2.5) 



where the disorder term comes from the dimensional analysis of the last term in eq. (2.1), which scale like IS}/ 2 N/R d / 2 . 
The Flory free energy for a Gaussian chain is then simply given by 

R 2 Nb 2 1/2 N 

^quench - ^2 + #T - A Rdfr ' y 2 - b > 

The first two terms corresponds to the stretching and compression entropy of the Gaussian chain. The third term 
comes from the quenched disorder. It is interesting for the discussion below, that the balance of the stretching 
term and the disorder term yields a critical value for the disorder A c cx N~ 2+d / 2 , whereas the comparison of the 
compression term and the disorder term provides the localized radius R ~ b(b d /A) 1 /*- 4- ^, which agrees with various 
calculations of the different type JjJ [l9| . This result for Gaussian chains states simply that the size of the chain is 
entirely determined by the disorder. The chain size does no longer depend on the chain length N, i.e., it is "localized" 
in a typical volume A. 



The corresponding Flory free energy for excluded volume chains has the form 35 



R 2 Nb 2 N 2 .1/9 N 
Quench = -ln^ typ « _ + _+«_- A 1 / 2 -^ , (2.7) 

where the third term corresponds to the excluded volume interaction. Because for the long chain N a / 2 <C N a , where 
a = 2 — vd, the excluded volume term dominates and the disorder leaves the Flory exponent unchanged. 

What can be expected from dynamic time scales? This can be estimated by similar arguments. Now we obtain for 
our random - field model a straightforward estimation of the diffusion coefficient D. The characteristic confinement 
time for the chain in a trap 

JV-l 



t* tt TR exp{^2v{IL(s,t)}} , (2. 



s=0 



where tr is the maximal Rouse time. The confinement time t* has a broad distribution according to fluctuations of 
barriers between traps. After averaging over these fluctuations we have 



. N—l N-l 



T«exp{-]T ]T 5(R(M)-R( S ',*))} 

s=0 s'=0 

»{^—} , (2.9) 



w T fl exp<j 

where the second line is a direct a la Flory estimation of the S - function terms, provided that the size of of the chain 
R - bN u . 

By taking into account eq.(|]^) the estimation for D reads 

R 2 r A 

D w - 



D R exp{-^N 2 -" d } , (2.10) 
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— N 2 ^ d » 1 



(2.11) 



which goes back to the Machta's result jjj], O. Since v = 3/(d + 2) the "specific - heat" exponent a — 2 — vd 
(4 — d)/(d+ 2) and it can be seen that at d < 4 and reasonably large 2V 

A 

so that the chain center of mass is effectively pinned down or frozen. 

The inequality (2.11) can be found independently as a paraphrase of the well - known Harris criterion p8[ which ex- 
press conditions when the disorder completely dominates over the chain entropy. We have relegated this consideration 
based on the n - component field theory formulation |34| to the Appendix A. 

From the dynamical standpoint the chain's center of mass freezing leads to the ergodicity breaking. It is a question 
of first importance to formulate a more general dynamical approach which enable to study the slow dynamics and 
ergodicity breaking not only for the center of mass but also for higher Rouse modes. We will devote the rest of this 
paper to the developing of such approach and show how in the long time limit the ergodicity breaking for different 
Rouse modes appears. 



III. EQUATION OF MOTION FOR THE TIME CORRELATION FUNCTION 

A. Langevin dynamics 



In this section we give a general consideration of the Langevin dynamics of a polymer chain in the quenched random 
field. The dynamics of the chain is described by the following Langevin equation 



Co j t Rj(s,t) - eA s R {s,t) 



-H int {R(s,t)} 



SRj(s,t) 



5Rj{s,t)~ 

v{n(s,t)} = fj(s,t) , 



(3.1) 



where j labels Cartesian components, Co is a bare friction coefficient and the second order finite difference A s Rj(s, t) = 
R 3 {s + l,t) + Rj(s - l,t) - 2Rj(s,t). 

The Langevin problem in question is getting much more convenient for the theoretical investigation if we change to 
the MSR - generating functional representation |36| . The generating functional (GF) of our problem can be written 
as 



£{•■•} = / DRj(s,t)DRj(s,t) 

x exp|A intra R(s,i),R(s,i) 



R(M),R.(M)]} 



(3.2) 



where the intra - chain action is given by 



N-l 



A h 



R(s,i),R(s,t) = dtliRj(s,t) 



s=0 



Co—R j (s,t)-eA s R j {s,t) 



5Rj(s,t) 

and the action related with the quenched random field reads 



H int {R 3 (s,t)} + k B TCo iRj(s,t) 



N-l 



R(s,t),R(s,t)\ = X! / dtiR j( s ^) 



s=0 



6Rj(s,t) 



V{R(s,t)} . 



(3.3) 



(3.4) 



The expressions ( p.2[ ) - (3.4) correspond to a given realization of the random field V{R(s,t)}. Now we perform the 
averaging over all configurations of X^{R(s, t)} taking into account its Gaussian statistics (see eq.( |2.3| )). The resulting 
GF takes the following form 



(Z{---}) v = I DR 3 {s,t)DR J {s 1 i) exp { A intra R(s,t),R(a,t) 

d d k 



N-l N-l 



/f d k 
dtdt' / -— 1 k j kie^p{ik[R(s,t) - R(s',t')}}iRj(s,t)iRi(s',t') 



(3.5) 
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It can be seen from eq.([3j|) that the averaging over the disorder leads to the non - Markovian (i.e. non - local in time) 
renormalization of the friction coefficient (which is coupled with iRj(a,t)iRi(a' ,t')). This causes actually dynamical 
slowing down and ergodicity breaking which we will discuss below. 



B. Self - consistent Hartree approximation 



In order to handle the functional integral (3.5) we use the Hartree approximation. In this approximation the full 
MSR - action is replaced by the Gaussian one in such a way that all terms which include more than two fields R(s, t) 
and/or R(s,t) are written in all possible ways as products of pairs of R(s,t) and/or R(s,t) coupled to the self - 
consistent averages of the remaining fields. On the other hand in ref.p^] it was shown that the Hartree approximation 
is equivalent to taking into account Gaussian fluctuations around the saddle - point solution. The resulting Hartree 
action is a Gaussian functional with coefficients which could be represented in terms of correlation and response 
functions. The calculation of these coefficients is straightforward and details can be found in the Appendix B of 



ref.! 



The second and third virial terms in Ai ntr8 



as well as the term which is responsible for the 



R(s,t),R(s,i) 

non - Markovian renormalization of the friction coefficient are treated in the same manner as in the ref. j39f . After 
collection of all these terms the final Hartree GF reads. 



N-1N-1 r oc 

EE 

s=0 s'=Q 
N-l N-l 

EE 

s=0 s'=0 

N-l N-l r00 



dt 



dt 



dt' iRj(s, t)Rj(s , i )A(s, a ; t, i ) 
dt' iR 3 (s, t)Rj (s, t)\(s, a'; t, t') 



I Ay ^ 11 * rOO poo 

oEE/ dt dt/iRj(8,t)iRj^,l/)x{.8,^]tM 

s=Q S '=0 J -° C J -°° 



where 



X(s,s';t,t') = ^G( S ,s';t,t') J 7^fc 4 F(k; a, s';t,t') 



J (^Vk)F st (k; S , S ') 
d d kd d q 



N 

E 

s" = l 



(27T) 



2d 



k 2 w(k, q)F st {q; a', s")F st {k; s, s') 



and 



X (s,s';t,t') = A J ^k 2 F st (k; S: s'). 
In eqs.(^.7[) - (3.8) the response function 

G(s,s';t,t') = ^R(s',i')R(M)) 

and the chain density correlator 

F{k; s, a'; t, t') = exp *, o} 



with 



Q(a,a';t,t') = (R(s, t)R(s, t)) - (R(s,t)R(s',t')> 



(3.6) 



(3.7) 



(3.8) 



(3.9) 



(3.10) 



(3.11) 



whereas F st (k; s, a') stands for the static limit of (3.1C). The pointed brackets denote the self - consistent averaging 
with the Hartree GF which is given by eq. (p\6j) . 
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In general one should consider fluctuation dissipation theorem (FDT) violation which is well known in the context 
of glass transition phenomenon pZif . In our present consideration we are mainly interested in the freezing conditions 
as well as the anomalous diffusion at the relatively short times. This enables us to assume that the FDT and the time 
translational invariance (TTI) are valid, then 



G(s,s'- 7 t-t') = (keT)- 1 — Q{s,s';t-t') at t > t' 



(3.12) 



By employing eq.( 3.12| ) in eqs.( |3.6[ ) - (3.8) and after integration by parts with respect to time argument t', we obtain 
the following Hartree GF: 



(Z{"'))v = DK-D&cxpi E / dt dt' iR 3 (s,t)[CoS{t-t') + 9{t-t')T(s,s';t,t')}—R :j (s',t') 

dt' iRjfat) Cl(s,s') Rj(s' ,t) 



E 

s,s'=0 



dt 



N-l r oo 



dt / dt 



' iRj(s, t) [£ 5(t - t') + 6{t - t')T(s, s'\ t, t')} iR 3 {s', i')} 



E 

where the memory function 

T(s,s';t,t') = A J ^Lk 2 F(k; S ,s';t,t') 
and the effective elastic susceptibility 

Q(s,s') = e<WA s - J ^fc 2 («(k)-A) 



(3.13) 



(3.14) 



JV-l 



F st (k;s,s')-S ss , F st (k;s,s") 



s"=0 



N-l 



F st (k; S ,s')F st (q; S ", S ')~5 ss , ^ F st (k; a, s'")F st (q; s'" , s") 



(3.15) 



In eqs.( 3.13 ) - ( 3.15 ) we use for simplicity the un its w here fc^T = 1, so that the disorder parameter A has the 
dimensionality of volume. The memory function fl3.14| ) is responsible for the non - Markovian renormalization of 
the Stokes fric tion coefficient Co which arises from interaction with the quenched field V^(k). The effective elastic 
susceptibility (3.15) takes into account all non - dissipative (reactive) forces in the system: local spring interaction, 



renormalization of the second virial coefficient due to the random field V(k) as well as third virial term. 



C. Equation of motion 

The equation of motion for the correlation function 

C(s,s';t,t') = <R(s,i)R(s',0) 



(3.16) 



can readily be obtained from GF (3.13). The result at t > t' reads 

d 



&—C(s,s';t,t) - J2 n(s,rn;t)C(m,s';t,t') 



rn—1 
N-l r t 



+ E / T (s,m;t,T)^-C(m,s';T,t')dT = 



(3.17) 
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It is convenient to make the Rouse transformation [ftl[ 



N-l 



and 



2jt 



(3.18) 



C(s,t)=J2C(p,t)eM-isp) , (3.19) 

p=0 

where p = 2nj /N (j — 0, 1, . . . N — 1), i.e. we have used for simplicity the cyclic boundary conditions. After this 
transformation the eq.(3.17) is simplified and takes the form 



d 



d 



Cog- t C(p;t)+N j T(p,t-f) — C(p;t')dt' + n(p,t)C(p;t) = 



where 



and 



NT(p,t) 



Of 



+ 2 



2d+l n d/2 



N-l 

E 



cos(ps) 



n=a [Q(n,t)] 



N-l 



1 — cos(pn) 



"(p) = i^ 1 - cos p) ( v - A ) S\ + t \* E 



1 — cos(pn) 



^W^S r^o [Qst(n)]^ [Q s tW] 4 



In eqs. ( 3.21 ) and ( 3.22 ) the time dependent mean - square distance 

Q{s,t) 



i([R(M)-R(0,0)] 2 



2tt 



J2iCst(p)-cos(ps)C(p,t)} 

p=0 



as well as its static limit 



(3.20) 



(3.21) 



(3.22) 



(3.23) 



Qst(s) = [1 - cos(ps)] C st (p) 

p=0 



(3.24) 



make the whole equation of motion for C(p, t) self - consistently closed. In the course of derivation of eqs. ( |3.20| ) - ( 3.22] ) 
we have took into account that the segment - segment interaction is short - ranged, i.e. w(k) « v and u>(k, q) « w; 
we have used also the Rouse transformation of the chain density correlator, i.e. 



1 JV_1 f 
F(k;p; t) = — ^ cos(pn) exp j 

n=0 ^ 



-Q(n,t) 



The static limit (i.e. t — > 0) is evident from eq. ( |3.20| ) provided that the initial condition 



Co[g- t C( P ;t) 



C o G( P ,t^0+) 



d_ 

'N 



is taken into account. Then the static equation becomes 



[^(pF 1 = i(l-ooBp)-(»-A) -^E 1 ^^ 



^-1 N-n-1 

IV Jd+T7Z)d E E 



1 — cos(pn) 



(3.25) 



(3.26) 



(3.27) 



„=0 m=0 [Qst(n)] 2 [Q st («)] 2 

It is of interest that this equation is similar, with an accuracy of prefactors and shifting v — > v — A, to the variational 
equation, which we have derived in ref. In the static limit this shifting is an only consequence of the random field 
V(r) effect. 
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IV. DYNAMIC BEHAVIOR OF THE CHAIN 



We are no w in position to launch a more elaborate investigation of the chain dynamic behavior which is based on 
eqs.( 3.20| ) - (3.23). There are at least two subjects which can be studied: i) the anomalous diffusion on the interval 
between a microscopic characteristic time Td (see below) and the longest internal relaxation time [jl], ^2| ; ii) Rouse 
modes dynamical freezing at t — > oo. 



A. Anomalous diffusion 

The presence of the quenched random field restricts the motion of the chain already at the time interval, 

r d < t < t N 1+2 " (4.1) 

where Td is a crossover time when the disorder starts to show up (the value of Td will be discussed below) and ToN 1+2v 
is the maximal Rouse time with the Flory exponent v. This restriction manifest itself through the subdiffusional 
regimes (anomalous diffusion) which have been seen first by Monte Carlo (MC) simulation ||. 
Let us start from the general solution of eq.( 3.20 ). For the Laplace correlator 



/•OO 

C(p,z)= / dtC{p, t) exp(-zt) 
Jo 



this solution reads 



C(p,z) 



( +Nr(p,z) 



z + 



(4.2) 



(4.3) 



The calc ulat ion of T(p, z) is based on eq.(3.21), where the time dependent mean - square distance Q(s, t) at the time 
interval (4.1) is approximated by 



20 

f)!>.n : : ir ( - ) -g,,(.v) 



(4.4) 



where 9 — v j(\ + 2v) and Q s t(s) — b 2 s 2w . This form can be just ified by im plem enting simple scaling arguments for a 
pure (i.e. without disorder) model |39|j . The substitution of eq.(L4) in eq.( |3.21 ) leads to the following result: 



NT(p -> 0, t) = A 



A \ /tq\/5 
b d+2 ) V t 



where the constant 



A 



^ +2 r(^)r(f-^ + i) 
2*+2 n i v f (f + i) 



(4.5) 



(4.6) 



with r (x) stands for the gamma function and the exponent 



P 



d + 2 



2v + l 



< 1 



(4.7) 



where a = 2 — vd is the "specific - he at" exponent. 
The Laplace transformation of (L5) at tqz <g; 1 reads 



NT{p -+0,t)=A 



[ b d+2 J T 



1-/3 



(4.8) 
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In the case when the memory term is in excess of the bare friction coefficient, i.e. at t > ra, we can use eq.(l.S) in 
cq. ( |4.3| ) which after inverse Laplace transformation can be put in the form (see the ref. 1 36 ) 



C( P ,t) = C st ( P )J2 



fe=0 



Center of mass mean square displacement 



AA 



(A) 



r(fc/? + 1) 



(4.9) 



.(*) = i([R c . m .(<)-Rc. m .(0)] 2 
= lim -C(p,t)} 



(4.10) 



The substitution of eq.(LG) in cq.( 4.10| ) results in the leading term of the anomalous diffusion, i.e. 



where 



b b + 2 



(4.11) 



(4.12) 



In the course of derivation of eq.( 4.11 ) we have used the static equation ( 3.27 ), i.e. C s t{p)£l(p) = l/N 



It is easy now to estimate the crossover time after which the disorder starts to effect the diffusion (see eq.(4.1)). 
The condition for that, Co = J"J" d dt T(p — > 0, t) can be recast in a form 



b d+2 Co 



A At, 



(4.13) 



One can see that the anomalous diffusion exponent f3 does not depend from the strength of disorder, whereas the 
prefactor Vq decreases with increasing A. For a chain in the good solvent v — 3/(d + 2) and at d — 3 the exponent 
Psasn = 0.9. For a Gaussian chain v = 1/2 and /fcauss = 0.75, i.e. the subdiffusional exponent has the same value as 
in a polymer melt |3q| . Finally in the case of the globule state v — 1/3 and /^Globule = 0.4, i.e. the globule anomalous 
diffusion is suppressed down by the disorder at most. 



B. Center of mass freezing 



Let us consider now the large time center of mass diffusion. In this case the characteristic time interval 

t > t N 1+2u (4.14) 

and internal Rouse modes are alr eady relaxed. For this time regime the reasonable approximation for Q(n, t) has the 
following form (compare with eq.(|4.4|)): 



Q(s,t) = dDt + Q st (s) 



(4.15) 



where D is the full (not bare) diffusion coefficient, which is renormalized by the effect of disorder and should be find 
self - consistently. On the other hand the equation for the zero - mode diffusion coefficient has the form pty, B8, f43| 



D 



1 



iv[Co + iv/ 00 ^r(p = o,t)] 



(4.16) 



(we recall that in our units of measurement kgT — 1). Eq. ( 4.16| ) enables to find D self - consistently. By making 
use eqs.(3.21) and ( 4.15| ) in eq.( 4.16| ) we obtain the following result for the center of mass diffusion coefficient 



D = D R (l-AFx 



(4.17) 
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where Z?r = (Co A) 1 is the Rouse diffusion coefficient and 



-7~N — 



d d / 2 

2 d n d / 2 



N-l 

N E 



IQst(s)} 



d/2 



(4.18) 



It can be seen that at AJ-^ > 1 the center of mass diffusion is frozen and the system becomes nonergodic. The relevance 
of this result is twofold. First this is a particular case of the so called A - type dynamical phase transition which 
has been extensively discussed [Q in the context of the mode coupling theory. On the other hand if we substitute 
Q st (s) in eq. fl4.18| ) with its most representative term Q st w b 2 N 2 " we will find D « D R [l - cons^A/b^A 2 " 1 ^] . As 
a result we turn back to the Machta's formula or, m ore exactly, to its expansion up to the first order with resp ect to 
(A/b 2 )N 2 ~ l,d . This actually means that eq. ( 4.17 ) overestimate the freezing and one should rather treat eq. ( 4.17 ) 
as a crossover criterion for the weak ergodicity breaking transition in the ref. J45| sense. 



C. Rouse modes freezing and a two mode toy model 



Now we study the freezing or the ergodicity breaking of Rouse modes with p ^ 0. This phenomenon mathematically 
manifests itself as a bifurcation with respect to the non - ergodicity function which, in its turn, is a long time limit 
of the corresponding correlator |3^, Q. Let us define such persistent part of the normalized correlator (i.e. non - 
ergodicity function) as the long time limit 



f(p) = lim 



C(p,t) 



*->oo C st (p) 

The equation for f(p) can be easily obtained by taking the limit t — > oo in eq.( |3.2C ). The result reads 

M 



2d+ ^i 2L, [L(s)]i+1 



1 - f(p) 



where 



271 



L(s)= C st (p)[l-cos(q S )f(q)} 

p=2ir/N 



(4.19) 



(4.20) 



(4.21) 



The eq.(4.20) is a self - consistent equation for the non - ergodicity function f(p). In the v icinity of the bifurcation 
point the non - ergodicity function f(p) is small and we can expand the r.h.s. of eq. (4.2C) with respect to f{p). It 
is shown in the Appendix B that because of orthogonality the zero - order term of this expansion vanishes and we 
arrive at the so - called F\ 2 - model according Gotze's nomenclature p3. The extensive numerical analysis of the full 



cq.( 4.20| ) which is given in the next section reveals that the bifurcation of f(p) is continuous or of A - type. 

T o ga in a better insight into the Rouse modes freezing mechanism let us consider first asimplified version of 
cqs.(B4). This is a toy model which is based on the truncation of the full hierarchical eqs.(B4) on the level of two 



longest modes, j — 1 and j = 2. In this case the asymptotic form 7VC s t ~ p 



-l-2v 



in order to calculate the coefficients in eqs.(B4). As a result the toy model equations for /(l) 
be recast in the following form 



where p w 2-irj/N <g; 1 can be used 



/ and /(2) = g can 



where the coefficients 



1-/ 
9 

1-5 



j2~vd 



Ai/ + ei/.g 
A 2 g + e 2 f 2 



(4.22) 



Ai = AA^ 

AN 2 ~ ud 



AN 



2-ud 



ei = 



A 2 = 



22+iv 



(2 



_ 2l+2^^^y2 — vd 



(4.23) 



It should be mentioned that the important f 2 - term in the second eq. (4.22) comes from the truncation of the last 
sum in r.h.s. of eq.(B4). 
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It is readily seen that in the vicinity of the critical point, = 1, the coefficient Ai = 1 + a, where a <C 1 and mode 
amplitudes have the following forms: / « crf + and g « <J 2 g+, where / + and g + are some constants. The substitution 



of these forms in eq.(4.22) leads to the solution 

fir) = ° 

g(°) = ° 2 y^ 2 ' (4 - 24) 

where it is important that A| < 1. 

As a result the trivial solution, / = g = 0, bifurcates at the critical point Ai = 1, so that the f-mode goes linearly 
and the g-mode - quadratically with respect to a. It is obvious that close to the critical point (i.e. a <C 1) there is no 
effect of the g-mode on the f - mode. On the other side g - mode bifurcates only as a result of f - mode bifurcation. 
In this respect one can say that the Rouse mode freezing follows the "host - slave" scenario. In Sec.V we will show 
that this scenario holds true for the whole numerical solution. 



NUMERICAL ANALYSIS 



In this section we present the numerical solution of cqs.(4.2C) - (4.21) in the full range of the Rouse mode index j 
values and for increasing values of the disorder strength A which here acts as a control parameter. As it usually is in 
the mode coupling theory fl33[ the full information about the static correlator, C s t{p), is a necessary prerequisite for 
the non - ergodicity equa tion study. In this respect, considering a chain of given length, we have numerically solved 
the static equation ( 3.27 ) for C st (p), where the virial coefficients and the disorder strength A are given. By making 
use of the Fast Fourier Algorithm we have implemented the bisection procedure between two trial profiles of C st (p) 
until the convergence to the final solution is achieved. This method has been recently used in the different context 
p2j where it has enabled t o con sid er ch ains of length up to N — 2 8 . After that we use C s t(p) as a static input for the 
non - ergodicity equation ( 4.20 ) - ( 4.21 ). This equation is solved for the chain length N = 128 in much the same way 
as it is described above for the static calculation. 



A. Bifurcation diagram 



We have found that for the small values of the disorder strength A the only solution of eq. ( 4.2CQ turn out to be the 



trivial one, i.e. f(p) = 0. As the disorder strength increases above a critical value A c , we observe that the first and all 
other modes become frozen simultaneously, i.e. they are characterized by a non vanishing value of the non-ergodicity 
function f(p) at the same A c . The resulting bifurcation phase diagram is shown in Fig.|l|. As may be seen from Fig.|l| 
all modes bifurcate continuously (A-type) , but bifurcations of higher modes (J = 2, 3, . . . ) are more smooth compare 
to the first mode bifurcation. This is qualitatively consistent with the result of the toy model analysis from the Sec. IV 
C. Moreover, one can see that the higher the Rouse mode index the more smooth is the bifurcation. This creates 
some numerical difficulties in precise location of the higher modes critical point A c . As soon as the accuracy of a 
non - vanishing value of f(p) (or the resolution) is not high enough the bifurcation diagram looks as if there were 
a subsequent mode freezing (see Fig.^]). We will show in the Sec.V B by the careful analysis of the finite resolution 
problem that this is an illusory effect and all modes do freeze at the same critical point. 

The critical value A c w 0.13 should be correlated with the radius - of -gyration diagram in Fig||. It can be clearly 
seen that at disorder strength comparable with this value the system is approaching the radius of gyratio n wh ich 



corresponds to the globular phase |42j . We recall that the static input information which is embraced by eq. ( 4.2C ) is 
determined by the effective virial coefficient v c g = v — A, where v is a bare second virial coefficient. That is why the 
bifurcation diagram in Fig.|l| corresponds to the Rouse mode freezing in the globule phase. 

It is interesting to elucidate how the Rouse modes freeze in the coil state. For this end we ha ve d riven the system 
to the coil state by increasing the value of the bare virial coefficient up to v — 0.5 while the eq.( [l.20] ) is solving. The 
result of the mode freezing is shown in Fig.^ and is compared with the previous case (where v = 0). It can be seen 
that the freezing of the modes in the coil state at least for small mode indices (0 < j < 10) occurs at higher values 
of A. This is intuitively clear and is consistent with our finding in Sec. IV A that the anomalous diffusion in the coil 
state is less affected by the disorder. 
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FIG. 1: Bifurc ation diagram for the non - ergodicity function according to 
equation (4.20). The calculation refers to a chain of N = 128 monomers and 
the bare second virial coefficient v = 0. We show only the smaller mode index 
values, i.e. j — 1, .., 5. The freezing of the modes appears above a critical 
value of A cr w 0.13. 




FIG. 2: 3D - Bifurcation diagram. The illusory Rouse mode successive 
freezing is a result of a finite resolution. In this case the accuracy of a non - 
zero value of f(p) or the resolution h — 10 -4 . 



B. Finite resolution and chain length study 



We have mentioned in the Sec.V A that because of very smooth bifurcation of higher modes it is numerically not 
trivial to locate there critical point A c . This location becomes very sensitive to the accuracy with which we measure 
a first non - vanishing value of f(p) while it bifurcates. We call this resolution h, and the FigJ| shows that as the 
chain length N increases and the resolution h is getting more fine (i.e. h decreases) the measurable critical points of 
the different modes merge each other. This proves that theoretically all modes freeze at the same critical point A c . 
Practically, since any experiment has a finite resolution one can presumably see that modes freeze sequentially (see 
Fig-D- 
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Disorder Strength, A 

FIG. 3: Radius of gyration as a function of disorder 
strength A for the chain length N = 128 and v = 0. 




Disorder Strength, A 

FIG. 4: Comparison of freezing diagram for coil and globule 
states. Two curves correspond to two different conditions with 
different bare virial coefficients, namely v = (solid line) and 
v = 0.5 (dashed line). 



Finally, in order to study how the critical disorder parameter scales with a chain length we have plotted (see Fig.g) 
the value of A c for different N. This numerical study, which the Figj6| depicts, reveals that the scaling law has the 
form: A c = o/V -7 , where the prefactor c increases with the solvent quality but 7 « 0.25 and is universal. This shows 
that A c for non - zero Rouse modes scales differently than for the zero - mode, which follows the Harris freezing 
criterion (see the Sec. IV B). Such scaling law obviously indicates the importance of the Rouse modes coupling. 
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FIG. 5: A c as a function of chain length N and the resolution h: a) h = 10 3 , b) h = 6.5 x 10 5 . 
Symbols correspond to the different modes: (■) j = 1; (*) j = 2; (x) j — 3; (+) j = 4 
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1000 10000 
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FIG. 6: Finite size analysis shows how the critical A c for the non - zero 
Rouse modes depends from the chain length N. Symbols correspond to 
the different solvent quality: (+) v = 0; (x) v = 0.5. This reveals 
the scaling law A c oc iV~ 7 , where 7 » 0.25. 



VI. CONCLUSION 



We have studied the dynamics of the self - interacting polymer chain (with second and third virial terms) which also 
experiences a quenched Gaussian random field. Consideration is based on the MSR - generating functional method 
and self - consistent Hartree approximation which allow to derive the equation of motion for the bead-bead time - 
dependent correlation function C(s, s'; t, t'). For the time interval less then the maximal Rouse time, t# <~ tqN 1+2 ", we 
have found an anomalous diffusion regime and have calculated explicitly the correspondent subdiffusional exponent. 
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At the larger time interval, t > tr, the ergodicity of the system is getting broken upon the disorder strength A 
increases (i.e. the chain becomes frozen). We have derived the corresponding equation for the non - ergodicity 
function f(p) and have solved it numerically. This solution shows that with increasing A all modes freeze at the 
same critical point, even though the higher modes bifurcate much more smooth. The center of mass (or zero - mode) 
dynamical freezing is governed by the Harris criterion, (A/b d )N 2 ~' /d > 1, whereas the critical point for the non - zero 
modes freezing is scaled as A c oc TV -7 . As it is found numerically 7 « 0.25 and even does not depend on the solvent 
quality. This universality is presumably the result of strong mode coupling effects. On the other hand it is of interest 
that the freezing line for coil and globule states are different with the coil becomes frozen at smaller disorder strength 
A (see Fig.g). 

It should be emphasized that the whole consideration in this paper is based on the validity of FDT and time transla- 
tional invariance (TTI). It is sufficient for the discussion of Rouse modes freezing close to the critical point A c . For the 
investigation of the gyration radius freezing well below A c it is necessary to discuss "R? g (t) — (1/A) Y^=o < R 2 { s i t) > 
in the aging regime where FDT and TTI do not hold poj] . We will return to this study in following publications. 
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APPENDIX A: HARRIS CRITERION 



It was argued presumably for the first time by Machta |l2j, |l3| that even though the disorder has no effect on the 
Flory exponent it can influence the whole spatial distribution of the chain at d < 4. Before him some authors [|l6| [lTj 
believe that disorder is completely irrelevant for the chain configurations provided that v > A. In this Appendix we 
paraphrase the Harris criterion argumentation given in ref. [ p8[ and apply it to a self - avoiding chain problem. 

Let us start from the field theory representation for a self - avoiding chain in the presence of a quenched random field 
V(r) with the Gaussian distribution and the second moment given by eq.(2.3). We will use the usual n - component 
tpa - field theory Hamiltonian |34fl which for our problem has the following form: 



I uz 1 1 

11 1 ddr \ Yd £ ( v ^w) 2 + r ^ v ^ (r) + r 

L a—1 a—1 



- 2 V(r) 



a=l 



(Al) 



where ip a ( r ) (a = l,2,...n) is a n - component field and r is the chemical potential conjugated to N. In this 
representation V(r) acts as a random fluctuations of the effective transition temperature: r — * t — V(r). 
The local minima configurations are given by the saddle - point equation 



V 2 VaM + [t - V(r)\ Mr) + vrp a {r) £ ^ 2 b (r) = 



(A2) 



6=1 



The Harris criterion should follow from the condition that in the average the quenched fluctuations of V"(r) dominate 
over the chain's entropy. Let us take a region of a linear size L larger then chain's size, i.e L 3> R. The typical value 
of the frozen fluctuations in this region reads 



V 



type 



A 1 / 2 

JJJ2 



Similarly to ref.[^4| let us assume that the field has 0(n) internal symmetry and can be parameterized as 

ipa{r) = n a 4>{r) , 



(A3) 



(A4) 



where ft is a unit vector in the space of the n - vector model, i.e. X]™=i n a = !■ The vector field ip a corresponds to 
the polymer density in the following way: p(r) — (1/2) J22=i ^a( r ) = (l/2)'/ ,2 ( r )- Then the field averaged over the 
whole L - region obeys the equation 



V 



type 



V(f> = 



(A5) 
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where (j) ~ p. If the fluctuations Vtypc are strong, i.e. Vt ypc 3> t, the by making use of eq.(A3) we get 

A i/d 



L < 



r 2/d 



On the other hand L 3> R ~ for where r ~ 1/iV (i.e.i? ~ iV"). Combining this with eq. (A6) yields 



(A6) 



(A7) 



The "specific - heat" exponent a = 2 — dv = (4 — d)/(d + 2) is positive at d < 4 and under condition ( A7) the disorder 
substantially affects the chain's statistics. This recovers the Harris criterion which we have discussed in Sec. II B. 



APPENDIX B: CONNECTION WITH F 12 - MODEL 



It is instructive to show how our general eq.(4.2C) is related with Gotze's F12 - model [33 . This reduction lays a 
solid mathematical basis for our analysis. 

First of all we will need the following orthogonality conditions which can be readily obtained by direct calculations: 



N-l 



cos(ps) = N5. t 



s=0 



iV-1 



cos(ps) cos(qs) 



pO 



-NS V 



(Bl) 



(B2) 



s=0 



N-l 



V cos(ps) cos(<zs) cos(ks) = — { S x K ' p+q | S , K ' q - p 1 



if q > p 
if q < p 



(B3) 



Now we can expand r.h.s. of eq.( 4.20| ) with respect to the non - ergodicity function keeping in mind that close to the 
bifurcation point f(p) is small. Then performing this expansion up to the second order and taking into consideration 
eqs.dBl]) - (pf) yields 



m 
1 - f(p) 



d? +1 (d + 2) N 2 C 2 st (p) 



i+2 



f(p) 



d + 4 



8c Bt (p) e; 



^2tc 



=2tt/N 



Cst(g) 



2tt 



2tt 



x [ C st (q)C st (p + q)f(q)f( P + q)+ £ C st (q)C st (\p - q\)f(q)f(\p - q\ 

q=2n/N q=2lt/N 



(B4) 



where E' ' q=2-n/N means that term p = g is dropped out. It is of interest that because of eq.(|Bl) the zero - order term 
in this expansion is vanished, so that eq.(B4) always has the trivial solution f(p) = 0. The cq.(Bl|) has the functional 
structure of F12 model which has been studied by Gotze j33| . 
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